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Abstract 



By means of the LSDA+U method and the Green function method, we 
investigate the electronic and magnetic properties of the new material of 
SrsCaRe3Cu4024. Our LSDA+U calculation shows that this system is an 
insulator with a net magnetic moment of 1.01 ji^/i.M.^ which is in good agree- 
ment with the experiment. Magnetic moments are mainly located at Cu 
atoms, and the magnetic moments of neighboring Cu sites align anti-parallel. 
It is the non-magnetic Re atoms that induce an orbital order of d electrons 
of Cu atoms, which is responsible for the strong exchange interaction and 
the high magnetic transition temperature. Based on the LSDA+U results, 
we introduce an effective model for the spin degrees of freedom, and inves- 
tigate the finite-temperature properties by the Green function method. The 
obtained results are consistent with the experimental results, indicating that 
the spin-alternating Heisenberg model is suitable for this compound. 
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Perovskite transition-metal oxides have been investigated intensively [1]. The discoveries 
of high-T c superconductivity [2] and colossal magnetoresistance (CMR) [3] in this family of 
materials pose significant challenge to theoretical understanding of the strongly correlated 
electron systems. 

In certain parameter regimes, strong Coulomb repulsions make electrons localized and the 
materials behave as Mott insulators. In this phase, low-energy excitations are dominated by 
the superexchange interaction [4], and the collective excitations of spin degrees of freedom 
govern the physics. Undoped high-T c superconductors (HTC) are such typical examples. 
Their properties are well described by a spin- 1/2 antiferromagnetic (AFM) Heisenberg model 
on a square lattice. It is widely believed that properties of the undoped systems are closely 
related to the origin of high-T c superconductivity [5] . 

It is also well known that the charge distribution in perovskite transition-metal oxides 
is influenced strongly by the crystal field, i.e. the degeneracy of the 3d electrons is lifted by 
lattice distortions of the perovskite structure in various ways. Therefore, the spin, charge 
and orbital degrees of freedom are intervened with each other, which make the physics of 
the perovskite transition-metal oxides very rich. These features are shared by materials in 
which CMR are observed [6]. 

As will be revealed in this Letter, the recently synthesized material Sr 8 CaRe 3 Cu4024, 
which forms the cubic perovskite structure as shown in Fig. 1(a), is another member of this 
family of materials [7]. The new material is similar to the undoped HTC in the following 
aspects: (1) perovskite structure, (2) magnetic insulator, (3) Cu surrounded by O. However, 
it is different in the points that this material is three-dimensional and has nonzero magne- 
tization in the ground state M ~ 0.95/xb per formula unit (f.u.) Also, this compound is 
peculiar since ferromagnetic (FM) cuprates are very rare, and has an unusually high mag- 
netic transition temperature T c ~ 440 K. Curie temperatures T c of known FM cuprates 
are usually quite low. For example, the T c of La4BaCu 2 Oio, K 2 CuF 4 and SeCu03 are 5, 
6.5, and 26 K, respectively [8]. Therefore, it is interesting to see why a strong FM state 
is realized in this material. Although the response of the compound to carrier doping is 
also very interesting, we will concentrate here on the undoped material. It is shown that 
both the electron-electron correlation and the coupling between orbital and spin degrees of 
freedom are important for this compound, which results in an spin-alternative Heisenberg 
Hamiltonian. 

We calculate the electronic and magnetic structures of Sr 8 CaRe 3 Cu4024 by using the 
WIEN2K package [9] , which is an implementation of the density-functional APW+lo method 
[10]. The MT sphere radii of 2.0, 2.0, 2.0, 1.9 and 1.5 a.u. are chosen for the Ca, Cu, Sr, Re, 
and O atoms, respectively, with i?X max =7.0, which results in about 5,200 LAPWs and local 
orbitals per cell. We use 700 k points in the Brillouin zone. As for the exchange-correlation 
potential, we adopt the standard generalized gradient approximation (GGA) [11]. Our 
GGA calculation predicts this material as a metal inconsistent with the experiment, which 
indicates that the electron-electron interactions are important for this compound. Thus, we 
use the LSDA+U method [12] to take them into account, and adopt £7=10 eV and J=1.20 
eV for the d orbital of Cu [13]. 

In order to investigate the magnetic structures and interactions of Sr 8 CaRe 3 Cu4024, we 
calculate both AFM and FM configurations for Cul and Cu2 moments. The results are 
shown in Table I. Magnetic moments are mainly located at Cul and Cu2 sites in both 



2 



configurations with the magnitude almost independent of the configuration. The magnetic 
moments at Ca, Sr, 01 and 03 are smaller than 0.001 /jl-q. Thus, they are negligible com- 
pared to the moments at Cu. 02 carries a small but nonvanishing moment due to large 
hybridization with Cul and Cu2. Our numerical results show that the AFM configuration is 
the ground state with net magnetic moment M=1.01 /ie/f-u. This is quantitatively consis- 
tent with the experimental one. Also, the ground state is an insulator with the energy gap of 
1.68 eV consistently with the experimental result [7]. It should be noted that the magnetic 
moment at Re atom is very small (less than 0.004 /i B ) in both configurations. This feature 
is contrasted with the ordered double perovskite A 2 FeM06, where due to the hybridization 
with d orbitals of Fe, the non-magnetic ion M possesses a large magnetic moment and plays 
a significant role in the magnetic properties [14,15]. We also performed calculations for U=8 
and 12 eV to cover a wide range of generally accepted values of correlations. Numerical 
results show that the magnetic moments are not sensitive to U . With decreasing £7, the en- 
ergy difference between AFM and FM (AE) increases slightly. For example, the calculated 
AE for U=8 eV is higher than that of U=W eV by about 20%. 

Figure 2(a) shows the calculated charge density plotted in (010) plane of Fig. 1(a). The 
charge density along the Re-03 bonds is larger than that along Cu2-03, indicating a stronger 
bond for Re-03 than that for Cu2-03. Therefore, the 03 atom is attracted toward Re. This 
induces a Jahn- Teller distortion in the oxygen octahedron centered at Cu2 consisting of 
four 03 and two 02 atoms with the bond length of Cu2-02 smaller than that of Cu2-03. 
Consequently, the e g orbital of Cu2 splits into d 3z 2_ r 2 and d x 2_ y 2. The former which points 
to 02 has a higher energy than the latter. Thus as shown in Fig. 3(b), the partially occupied 
orbital in Cu2 is the minority spin (spin-up) rf 3z 2_ r 2, and the charge density distribution of 
Cu2 is anisotropic as shown in Fig. 2(a). In contrast, three crystallographic directions of 
Cul are completely equivalent, and the oxygen octahedron centered at Cul site is free of 
any distortion. Therefore, the e g orbital of Cul is still fully degenerate, and all d orbitals 
except for the minority spin (spin-down) in the e g orbital are fully occupied as shown in Fig. 
3(a). In addition, the 2p states of 02 distribute mainly in the energy range from -7.0 to 0.0 
eV, and are almost full-filling. 

The spin density shown in Fig. 2(b) clearly indicates that the magnetic moments of 
the Cul and Cu2 are carried mostly by e g and d 3z 2_ r 2 orbitals, respectively. In Fig. 4, 
we schematically show the partially occupied d orbitals of Cul and Cu2 and the almost 
fully occupied p z orbital of 02 for the material. It is interesting to observe that an orbital 
ordering appears at Cu sites. Namely, the minority spin d 3z 2_ r 2 orbital of Cu2 and the 
e g orbital of Cul are pointing to the 02 as shown in Fig. 4. Both of them are less than 
half-filling. Meanwhile, the fully occupied p z orbital of 02 points to the neighboring Cul 
and Cu2 sites. Therefore, the spin-up and -down p z orbitals of 02 strongly overlap with 
the spin-up d 3z 2_ r 2 of Cu2 and the spin-down e g orbital of Cul, respectively, to form a 
rather strong pda hybridization. This results in a strong exchange interaction between the 
magnetic moments at Cul and Cu2 and a high transition temperature. 

In order to further clarify the effects of Re, we perform a calculation for an artificial 
structure, where the Jahn- Teller distortion in oxygen octahedron centered at Cu2 is removed. 
We find that the orbital ordering still survives, and the results are very similar to that of the 
real structure. We then change Re to Technetium, which is the element at the same column 
in the periodic table. Our calculation shows that the partially occupied orbital of Cu2 
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changes from d 3z 2_ r 2 to d x 2_ y 2. Consequently, the spin and the orbital are almost decoupled, 
and the spin interaction becomes much weaker. As a result, the energy difference between 
AFM and FM configurations is reduced greatly. Therefore, the Re atoms play an important 
role for the orbital ordering to cause the unusually high T c , although their magnetic moments 
are very small. 

Since the magnetic moments are almost localized at Cu sites, the effective model for 
the spin degrees of freedom of this material is expected to be a Heisenberg model, where 
neighboring localized spins at Cu sites interact antiferromagnetically. Although, the simple 
ionic model (Cu 2+ -0 2_ or Cu 3+ -0 2_ ) is not exact due to the large hybridization between Cu 
and O, we still take spins at Cu sites as a good quantum number similar to the case of the 
undoped HTC. We consider the following two cases: (l)5i=l, 5 2 =l/2 and (2) S r 1 =S r 2 =l/2, 
where Si denotes the spin length at Cm site. The magnetization in the unit cell for the 
AFM state and that of the FM state are obtained as (1) M AFM =l/2, M FM =5/2 and (2) 
M AF m=1, M fm =2. The LSDA+U result is M AF m=1.01 Mb ^0.5, M FM =5.01 /i B ^2.5 as 
shown in Table I. Hence, case (1) is reasonable. Thus, the effective Hamiltonian becomes 
the following: 

H = JeK^Si'^Si+z, (1) 

i P 

where J e ff is the effective exchange interaction, i runs over Cul sites, and p denotes the 
unit vectors of the unit cell in the lattice shown in Fig. 1(b) (p=±x, ±y, ±z). Here, Si and 
s i+ p denote the spin operator of S—l at Cul-site(z) and that of 5=1/2 at Cu2-site(i+|), 
respectively. 

Since this model is on a bipartite lattice and belongs to the family in which the Marshall- 
Lieb-Mattis theorem holds [16], the ground state of this model is proven to have spin 
S = l/2(=l/2x3-l) in any size of systems. This results in spontaneous magnetization M—l 
/iB in the bulk limit, and well explains the experimental result (M ~ 0.95/xb [7]). 

Since the difference of the magnitude of magnetic moments at Cu sites between FM 
and AFM states is small as shown in Table I, the energy difference between FM and AFM 
A£'=0.036 Ry=5,679 K can be assigned to the exchange interaction. By diagonalizing 
the Hamiltonian (1) in the unit cell, we evaluate the energy difference between FM and 
AFM configuration as 8J e ff- Hence, we estimate J e Q=AE/ 8=710 K. Using the above ef- 
fective Hamiltonian, we investigate magnetic properties of Sr 8 CaRe 3 Cu4024. In order to 
take quantum fluctuations into account, we apply the Green function method [17]. In this 
method, we construct sixteen Green functions of Cul-spins and Cu2-spins. We use Tyab- 
likov's decoupling [18] to obtain closed coupled equations for Green functions. Following 
Callen's scheme [19], we obtain (s z ) and (S z ) as a function of spin- wave dispersion relations 
Ek- Here, Ek is also a function of (s z ) and (S z ). Thus, we solve these coupled equations 
numerically to obtain (s z ) and (S z ) at finite temperatures. 

The result of the Green function method for the spontaneous magnetization is shown 
in Fig. 5. In the low-temperature limit, spontaneous magnetization M obtained by the 
Green function method is about 1/xb, which is consistent with the exact result M—l fi B 
and the experimental result M ~0.95 /ie- Magnetic transition temperature T c is calculated 
as T c = 448 K by the Green function method, which is close to the experimental result 
T c ~440 K. Thus, the magnetic properties known so far in Sr 8 CaRe 3 Cu4024 are explained 
by the spin-alternating Heisenberg model defined in eq.(l). 
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Here, we summarize the mechanism of ferrimagnetism and the high magnetic transition 
temperature T c for Sr 8 CaRe 3 Cu40 2 4. The mechanism of the AFM coupling is the superex- 
change [4]. The large exchange constant J e ff is due to the large overlap between Cu and 02 
orbitals, and to the localization of Cu-spins by strong correlations. This situation is similar 
to the case of the undoped HTC. In contrast to the undoped HTC, the number of sublattice 
sites is different. This causes ferrimagnetism. Imbalance of the spin length between Cul 
and Cu2 sites would also be a reason for stabilization of the ferrimagnetic state. 

In summary, using the LSDA+U method and the Green function method, we have inves- 
tigated the electronic and magnetic properties of Sr 8 CaRe 3 Cu4024. Our results show that 
magnetic moments are almost located at Cu sites, and the moments of Cul and Cu2 align 
anti-parallel due to the superexchange interaction [4]. The ground state is a ferrimagnetic 
insulator with a net magnetic moment of 1.01 /xs/f.u., which is in good agreement with the 
experimental results. We find that the non-magnetic transition metal Re plays an important 
role in determining the magnetic properties of this material. Namely, Re induces orbital or- 
derings of the d orbitals in Cu. This results in the strong exchange coupling between the 
magnetic moments of Cu and the unusually high T c . Based on the LSDA+U results, we have 
introduced an effective model for the spin degree of freedom, and investigated the magnetic 
properties by the Green function method. The ground-state magnetization of this model is 
shown to be exactly 1 /xb, which is consist with the experimental result. The obtained T c is 
also consistent with the experimental one. Therefore, the spin-alternating Heisenberg model 
is suitable for describing the magnetic properties of this compound. 
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and for valuable discussions. Dr. M. Arai and Dr. M. Katakiri are appreciated for technical 
helps. This work was partially supported by Japan Society for the Promotion of Science 
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TABLES 



TABLE I. Calculated total energy, Fi tot relative to the energy of AFM configuration in units of 
Ry, the total magnetic moments per unit cell fi tot , the magnetic moment at Cul, Cu2 and 02 in 
units of hb> 







P-tot 


Cul 


Cu2 


02 


FM 


0.036 


5.01 


1.15 


0.84 


0.14 


AFM 


0.0 


-1.01 


1.09 


-0.81 


0.07 
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FIGURES 

FIG. 1. Crystal structure of SrsCaRe3Cu4024. (a) Unit cell, (b) Unit cell of the effective 
model. 

FIG. 2. Contours for charge density in the (010) plane with interval 0.03 e/bohr 3 . (a) Total 
charge density, (b) Spin density (spin- up charge density minus spin-down charge density). The 
dotted lines are negative contours. 

FIG. 3. Projected density of state with Fermi energy at zero, (a) For minority spin (spin-down) 
d orbital of Cul. (b) For minority spin (spin-up) d orbital of Cu2. 

FIG. 4. Schematic picture for the spin and orbital orders of SrsCaRe3Cu4024 deduced from 
the LSDA+U calculation. The gray symbols denote partially occupied d orbital, and the dark ones 
denote full-occupied p orbital. The arrows denote the directions of magnetic moments on Cu sites. 

FIG. 5. Temperature dependence of the spontaneous magnetization. Solid curve denotes the 
Green function result. Squares and the dotted line are experimental results in Ref.[7]. 
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